Effects of electron-phonon coupling on Landau levels in graphene 
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We calculate the density of states (DOS) in graphene for electrons coupled to a phonon in an 
external magnetic field. We find that coupling to an Einstein mode of frequency uie not only shifts 
and broadens the Landau levels (LLs), but radically alters the DOS by introducing a new set of 
peaks at energies E„ ± uie, where E n is the energy of the nth LL. If one of these new peaks lies 
sufficiently close to a LL, it causes the LL to split in two; if the system contains an energy gap, a 
LL may be split in three. The new peaks occur outside the interval (— u>e,uie), leaving the LLs in 
that interval largely unaffected. If the chemical potential is greater than the phonon frequency, the 
zeroth LL lies outside the interval and can be split, eliminating its association with a single Dirac 
point. We find that coupling to an extended phonon distribution such as a Lorentzian or Debye 
spectrum does not qualitatively alter these results. 

PACS numbers: 73.22.Pr, 71.70.Di, 63.22.Rc 



I. INTRODUCTION 

Graphene, a single sheet of carbon atoms arranged in 
a honeycomb lattice, has aroused much interest in re- 
cent years for its many remarkable properties) 1 ^ Most of 
these properties follow from the fact that at low energies, 
the charge carriers obey the Dirac equation for massless 
fermions. Their dispersion relation is linear in momen- 
tum, with two bands, valence and conduction, that each 
form a cone called the Dirac cone, the tips of which meet 
at a point called the Dirac point. Angle- resolved pho- 
tocmission spectroscopy (ARPES), which provides a di- 
rect measure of the electronic dispersion curves, has ver- 
ified this band structure!^' From the dispersion, other 
interesting characteristics can be calculated: for exam- 
ple, an electronic density of states that increases lin- 
early out of the Dirac point, in striking contrast with 
the constant value found in conventional metals; and an 
optical conductivity that displays an intraband Drudc 
peak at zero frequency, followed by a sharp rise to a 
constant interband contribution at twice the chemical 
potential /zP^ These predictions, too, have been con- 
firmed, in scanning tu nnelin g spectroscopy (STS)^l!2land 
optical measurements QJIB] However, all of the above ex- 
periments have also demonstrated that the simple pic- 
ture of independent fermions must be modified by many- 
body renormaliz ations du e to electron-electronQSHEO anc j 
electron-phonorP l 6 * 9 * 10 ! 16 * ^ interactions. In particular, 
sharp kinks are observed in the linear dispersiorPHsl an j 
DOSp which have been identified as the result of coupling 
to a phonon. And optical measurement^^ find signifi- 
cant absorption in the region between the Drude peak 
and the constant background, a phenomenon thought 
to be at l east pa rtially caused by the electron-phonon 
interaction! 8 ! 18 ! 21 ^ Similar absorption absent from the 
independent-particle prediction^ has also been observed 
in bilayer systems 

These many-body effects may potentially be amplified, 



and new ones created, by applying a magnetic field to the 
system. In a field B perpendicular to the graphene sheet, 
the linear energy bands are replaced by a set of discrete 
Landau levels with spacing proportional to ^/|n|_B for in- 
tegral n. The effect of this q uantization on the DOS 24 
and optical conductivit y^ 25 * 26 ! have been calculated in an 
independent-particle pictur e, and those predictions have 
been verified in experiment j 1 l 1 ° l 2 '3 EIo| B u t we can expect 
these leading-order predictions to be modified by an as- 
sortment of many-body effects made possible by the more 
complex electronic structure of the quantized system. In 
particular, we can expect to observe a rich interplay be- 
tween the LLs and phonons, because for typical magnetic 
fields used in experiment, the LL spacing is on the scale 
of the phonon energies in graphene. Our goal here is to 
extend the independent-particle predictions by taking a 
first step, within a simple model, toward characterizing 
that interplay as it manifests in the electronic DOS. 

In a companion paper pD we showed that electron- 
phonon coupling induces a new set of peaks in the DOS, 
which we called phonon peaks, and that these peaks 
cause nearby Landau levels to split in two. We also 
showed that this can explain a disruption of the LLs 
at a characteristic frequency (which we identified with 
the phonon frequency) in STS dataP^ In the present pa- 
per, we complement those results by providing analytical 
calculations and considering a wider range of situations, 
focusing on descriptions of the effects of coupling rather 
than on quantitative predictions. We present the neces- 
sary formalism in Sec. [TTJ where we determine the loca- 
tions of all phonon peaks by calculating the self-energy 
in the case of coupling to a single Einstein mode of fre- 
quency loe- In Sec. |III| we show numerical results for the 
renormalized DOS in the case of zero chemical poten- 
tial, emphasizing the locations of the renormalized LLs 
and their interactions with the phonon peaks. We also 
present some results in the case of a Lorentzian distri- 
bution of phonons, showing that features are smeared 
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but qualitatively unchanged from those of the Einstein 
model. In Sec. IIVI we discuss the case of finite chemical 
potential, analyzing the locations of the phonon peaks, 
which become discontinous functions of B, and the pos- 
sibility of splitting the zeroth LL when [i is greater than 
uie- In Sec. |VJ we discuss two more scenarios in which 
electron-phonon coupling induces new structures in the 
DOS: the presence of an energy gap, and coupling to an 



acoustic phonon. We conclude in Sec. VI 



II. FORMALISM 
A. Governing equations and parameters 

In a magnetic field B, the bare electronic density of 
states with zero chemical potential consists of a line at 
each Landau levelP^l 



N baIe (cj) = ~N B9(W-\u\) 



S(oj + A) + 6(w - A) 



(1) 



where M n — sgn(n) yj2\n\B + A 2 is the energy of the nth 
level and W is a high-energy cutoff. The overall scale of 
the DOS is set by No = nh 2 v 2 and B, where the latter 
is in units of energy squared. To express B in units of 
teslas, we let B — > eBhvp (in SI units). We have allowed 
for the possibility of an energy gap A about the Dirac 
point, which splits the n = LL. Such a gap is closely 
related to an asymmetrybetween the electron densities 
on the two sublattices 1221221 Accounting for the cutoff W, 
if A = then the bare DOS has peaks at ui = M n for 
— Yb — n — Yb\ ' wri ere [^J is the largest integer 
smaller than x. If a finite gap opens, the zeroth peak is 
split into two, located at oj = ±A, while the other peaks 
are only shifted by the inclusion of A 2 in M n . 

In practice, the lines in the DOS are broadened into 
Lorentzians as S(x) — > ^ x 2+ r 2 ■ For simplicity, we will 
take r to be constant; in a more comprehensive calcula- 
tion, it would be determined by the imaginary part of the 
full electron self-energy, and it would depend on both ui 
and n. To orient the reader, we show the typical form of 
the bare DOS with finite T in the upper frame of Fig. [I] 
The bare DOS at finite B (shown in solid black) consists 
of a slowly varying piece, given by the DOS at B = 
(shown in dashed red), plus oscillations due to the mag- 
netic field. The slowly varying piece is contributed by 
the sum of tails of Lorentzians, which will be significant 
in what follows. 

As we shall see, the electron-phonon interaction rad- 
ically alters this set of peaks. The renormalized DOS 
is given by the convolution of the bare DOS with a 
Lorentzian broadened and shifted by the electron-phonon 
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FIG. 1. (Color online) The bare DOS with a magnetic field 
(solid black) and without (dashed red), along with the cor- 
responding real and imaginary parts of the self-energy cal- 
culated from those bare densities. Integers label the bare 
Landau levels. Here A = 50meV, oje = 200meV, F = 5meV, 
and the black curves are calculated for B = 15. 2T. Note that 
No has dimensions of inverse energy squared, so the quantity 
N/No has dimensions of energy rather than the density of 
states' usual dimensions of inverse energy. 



self-energy: 
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where E = Si +iE 2 is the self-energy One can write this 
simply as a sum of Lorentzians, 

N(u) = N Q B^-9(W -\u>- Ei(cj) + 

r-s a (w) 



[u - Ei (W) + fi - A] 2 + [r - E 2 (w)] 2 

r-s 2 H 
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(3) 



Alternatively, following Sharapov et al.p^one can eval- 
uate it to find 
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where Q = u> — Ei(w) + ^ + iT — iE2(w) and ip is the 
digamma function. This latter form of N is more conve- 
nient for numerical implementation. 

At zero temperature, the self-energy is given b} * 19 * 21 ! 



E( W ) 



1 

w 



r 

dvo?F(v) / 



e{-uj>) 



(5) 



where a 2 F is the electron-phonon spectral density. Note 
that, for simplicity, we will not take into account correc- 
tions due to the phonon self-energy as modified by inter- 
action with Dirac electrons, which have been previ ously 
calculated theoretically^^ and seen in experiment! 37 * 38 ! 

Equations Q and |5| are solved iteratively, obtaining 
the nth iteration JVWfrom EW, and E< n ) from TV^™" 1 ). 
We begin this process with the noninteracting DOS, , 
obtained by letting lo — » u + /xo in ATj, are , where is the 
bare chemical potential. This is equivalent to setting 
E = E<°) = and fj, = [i in Eq. ^ or Q. In subse- 
quent iterations, the noninteracting chemical potential is 
altered to the interacting one according to \i =_uo+Ei (0) , 
which preserves the total number of electrons 3 ^ and keeps 
fixed the value of the DOS at the Fermi energyPS Since 
we wish only to characterize the general behavior of the 
renormalized Landau levels and the locations of any new 
peaks, throughout most of this paper we stop after one 
iteration, which allows us to cleanly identify the effects 
of renormalization. We will discuss the effects of further 
iteration in Sec. MI Dl 

Again only to orient the reader, we show the typical 
form of the first iteration of the self-energy, calculated 
from the bare DOS, in the lower two frames of Fig. [I] As 
with the DOS, the self-energy consists of a slowly varying 
piece equivalent to the B = self-energy, plus oscillations 
due to the magnetic field. Here and throughout most of 
this paper, we adopt an Einstein model for the phonon 
distribution, given by a single phonon at frequency uje- 
a 2 F{v) — A8(v — oje), where A defines the strength of 
the coupling. As shown by Park et al.j^if is chosen 
to be 200meV, the energy of the E 2g mode in graphene, 
then in the absence of a magnetic field this model pro- 
vides a phonon-induced electron self-energy in excellent 
agreement with that from a full first-principles calcula- 
tion. Here we assume that it remains a good approxima- 
tion in the case of a finite field. Although the modes in 
the phonon spectrum may be weighted differently in the 
presence of a magnetic field — due to selection rules in the 
matrix elements that come from the vertices of the self- 
energy diagram — the Einstein model has already proven 
to yield qualitatively useful statements about the effects 
of electron-phonon coupling on LLsf^More importantly, 
since the self-energy simply sums over the couplings to all 
phonons, studying coupling to a single phonon provides 
a simple means of understanding and characterizing the 
effects of coupling to any distribution of them, regardless 
of their weighting. 



Although our goal is to analyze qualitative effects of 
coupling to a phonon, rather than to make quantitative 
predictions, we will for the most part use experimentally 
relevant values of our parameters. For the Fermi veloc- 
ity, we will always use the typical value vp = 10 6 m/s.^ 
In the Einstein model, we will almost exclusively use 
cue = 200meV, in accord with the results of Park et alP3 
To match to experiment J 3 -^ we used a coupling strength 
A ~ lOOmeV, but here we will generally use a larger 
value of 250meV in order to make qualitative effects 
more distinct. For the residual scattering parameter T, 
we will mostly use a value of 5meV, which lies in the 
middle of the rang e ~ .02-I0meV derived from trans- 
port measurements J22E1HH Ando has shown that a value 
of this ordei!^ is consistent with broadening due to im- 
purity scattering.^ For the high-energy cutoff we use 

where t is the nearest-neighbor hopping 
parameter; this ensures that the number of states in the 
Dirac cones equals the number in the first Brillouin zone. 
Specifically, we use W = 7eV, corresponding to the typ- 
ical value t ~ 3eV. However, the value of this cutoff is 
largely arbitrary. We can see from Eq. Q that it has only 
a marginal impact on the bare DOS, through the addi- 
tion of a constant term. And its impact on the dressed 
DOS is governed directly by the magnitude of the self- 
energy. The self-energy itself is more strongly affected 
by W, because W determines the limits of integration 
in Eq. (pH). Changing the limits of integration not only 
determines the total number of phonon peaks, but also 
the magnitude and slope of the slowly varying piece (i.e., 
the B — part) of the self-energy. The factor of 1/W 
only partially cancels this effect. However, both this and 
the overall effect on the dressed DOS can always be pre- 
cisely compensated for by changing A, since A appears as 
an overall prefactor in the self-energy. (We have verified 
this numerically.) Therefore, our model is independent 
of our choice of W in the sense that for a given value of 
W, any self-energy and dressed DOS can be arrived at. 
One must only keep in mind that the value of A is not 
meaningful in itself, since it depends on the choice of W. 



B. Analytical solutions and peak locations 

Before presenting numerical results for the dressed 
DOS, we provide an an analytical deduction of its peak 
structure, which will provide the foundation for the re- 
mainder of the paper. We begin by noting that the renor- 
malized Landau levels are located at the peaks of the 
Lorentzians in Eq. ([3]), where the term m~ Si (a;) +/x— M„ 
in their denominator vanishes. (Equivalently, the peaks 
lie at energies where the real part of the argument of ip 
in Eq. Q is a nonpositive integer.) The energies of these 
levels we denote by E n , which satisfies E n — T,i{E n )+n = 
M n . To determine the locations of any other peaks, we 
examine the structure of the self-energy. In the Einstein 
model, the imaginary part of E can be evaluated imme- 
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diately to find 

S 2 = 



ttA 



[N(uj — uje)Q{uj — <jJe) 



WN 

+N(lj + uje)0(-uj-uj e ) 



(6) 



Substituting N(u) = N^{uj) = N haxe (u + Ho) on the 

right hand side yields the first-order solution ■ To 
most easily see the peak structure, we let T = for the 
moment, such that after simplifying the Heaviside func- 
tions, £ 2 becomes 



^2 — 



— U(u -uj e + ho- A)0(- M o + A) 

+ 5(w -uj e + fx + A)0(-/i O - A) 
+ <5(w + w B + ho ~ A)6(hq - A) 
+ S(u> + lo e + no + A)6(po + A) 

+ 2J2 l 6 ^ -u E + Ho - M„)0(-Mo + M„) 



|n|<n m; 



S(lu + lu e + Ho - M n )6(p - M n ) 



(7) 



where n max 
by 



2D 



= — "P 
1 W 



Next, the real part of Y,^ 1 ' is given 
N^(oj') dJ 



No uj + uje — uj 1 
N^XcS) duo' 



A p [°°tl_^A ^ . ( 8 ) 

W J Nq uj — uj 'e — w' 



Again substituting N^°\uj) = Nb aTC (ui + ho) and letting 
r = 0, and taking careful note of when the supports of the 
delta functions in N lie within the range of integration, 
we arrive at 



= AB j 9(hq + A) | _ 
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where ni 



n=U2 + 



Mo 
2D 



l CJ + ^ - M n - UJ E 



(9) 
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Mo 
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TABLE I. Locations of peaks in E and N induced by the 
electron-phonon interaction, in addition to the logarithmic 
peaks at ui — ±uie- Each of these peaks lies at the position 
of a bare Landau level displaced by plus or minus the phonon 
energy. Top: the new peaks generated from the n = LL. 
Bottom: the new peaks generated from the LLs. 
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From these results we find that both £1(0;) and £2(w) 
have peaks at A — ho ± — A — /io ± wg, and M n — 
Ho i — in other words, the energy levels in the bare 
DOS, but shifted by plus or minus the phonon energy. 
We shall refer to these new peaks as phonon peaks and 
denote their energies as P m , where m is an integer. To 
distinguish this new set of integers from that labeling the 
Landau levels, we will add a subscript p, as in m p , when 
labeling energy levels. The values of P m are summarized 
in Table [TJ The upper frame shows the position of the 
m p — level, which depends on the relative values of Ho 
and A, a phenomenon which will be explored in Sec. |V A] 
The lower frame shows the positions of all other phonon 
peaks for the two cases ho < and Ho > 0. In each case 
the peak locations are divided into two sets, one contain- 
ing peaks lying to the left of —uje, one containing peaks 
lying to the right of +u>e — no peaks (including Pq) oc- 
cur inside the interval (— u>e, uje)- We shall find in what 
follows that this window, (— uje,^e), plays an important 
role in determining the general behavior of the DOS: in- 
side the window, the DOS consists simply of shifted LLs; 
outside of it, the DOS is radically altered by the appear- 
ance of the phonon peaks. In addition, the presence of 
the window affects the dependence of the peaks' position 
on the three parameters Ho, A, and \/2B: because no 
peaks can occur within the window, continuously vary- 
ing one of these three parameters so as to make a peak 
pass through the window's boundary at u>e or —uie forces 
the peak to jump discontinuously to the other side of the 
window. We will explore this phenomenon in detail in 
Sec. ITVAl 

Note that we cannot actually find N^{uj) from the 
self-energy just derived, because it would be a nonlin- 
ear functional of S functions. To self-consistently iter- 
ate the DOS and self-energy, we must begin with broad- 
ened S functions in Ab arc (o;). For S^, this has the ef- 
fect of simply broadening the S functions in Eq. Q into 

Lorentzians. For Ej 1 , we must reevaluate Eq. ([8]). Doing 



so yields 
E«(u,) = 



where 
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(w + w B + ^ - p) 2 + r 2 



(a; — lje + fJ-o — p) arccos 



V(A«o-p) 2 +r 2 
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(oj + uj e + - p) 2 + r 2 



Tin 



\u>—u>e\ 



V(p-Mo) 2 +r 2 

ue + po -p) 2 + r 2 ' 



(11) 



Each Lorentzian in N^°'{u>) has contributed a function, 
H p (uj), that depends on uj+uje and cj — u^, as we expect. 
In the limit r — > 0, the arccos terms in Eq. ( 10 ) reduce 
to Eq. (9), recovering the peak structure just described. 

However, for finite T, E^ contains logarithmic diver- 
gences at = ±ai£, which are not present in Eq. 
We can understand this by noting that the lo garithmic 
divergences are present in the B = case J^^XI and hence 
must come from the slowly increasing piece of the DOS, 
rather than from the peaks at the Landau levels. And 
since the slowly increasing piece arises from the sum of 
the tails of the Lorentzians, the logarithmic divergence 
exists only for finite F. Similarly, since N(lu) nowhere 
vanishes for finite r, E 2 always sharply changes from zero 
to the finite value — N(0) at w = ±u)e, as we can 
see from Eq. ffify. 

We are now in a position to understand the general 
form of E, as shown in the lower two frames of Fig. [T] 
for the simplest case of /i = A = 0. As discussed above, 
the self-energy for finite B (solid black curves) consists of 
a slowly varying piece (dashed red curves), constructed 
from the B — piece of the DOS (shown in the up- 
per frame), plus oscillations, constructed from the oscil- 
lations in the DOS. In fact, as we see in both the figure 
and Eq. (J6|, the oscillations in E 2 are precise images of 
those in the bare DOS, simply rescaled and translated 
by the phonon energy. Note the distinctive character of 
the peaks in Si and E 2 at ±lue- At those points, S 2 be- 
comes finite and picks up half an image of the n = bare 
Landau level, while in Ei the logarithmic singularity con- 
tributes to the P + peak picked up from the n = bare 
LL. Since no oscillations occur in the window (-wg, u>e), 
in that region the self-energy in the B > case cannot 
be distinguished from that in the B = case. 

Since the DOS has the form of oscillations atop a slowly 
varying function, each peak in the self-energy will intro- 
duce a peak in TV - ' 1 ' via that slowly varying piece. There- 
fore, from the above analysis we can conclude that the 
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FIG. 2. (Color online) Density of states and self-energy as 
a function of energy, for A — 50meV, loe = 200meV, B = 
15. 2T, and F = 5meV. The dashed red curve shows the bare 
DOS, while the solid black curve shows the DOS renormalized 
by E. Integers n label the dressed Landau levels. Integers m p 
label the phonon peaks, pointed to by blue arrows. 



renormalized DOS contains the following peaks: renor- 
malized Landau levels at energies E n labeled by integers 
n, which are solutions to E n — Ei (E n ) + p = M n ; phonon 
peaks at energies P m labeled with integers m p , which 
arise as images of the LLs in the bare DOS and always 
lie outside the interval (— oje,oje)', and phonon peaks at 
±Ub, which arise from a logarithmic divergence in Ei 
due to the abrupt opening of a new scattering channel 
at the phonon frequency. Depending on the values of 
Ho and A, the zeroth phonon peak Pq may or may not 
be split in two and may or may not coincide with the 
logarithmic divergences. Furthermore, as we shall see in 



Sec. IIIB each of the renormalized LLs may be split by 



the presence of the phonon peaks. 



III. DENSITY OF STATES FOR VANISHING 
CHEMICAL POTENTIAL 

A. Appearance of phonon peaks 

We begin our numerical analysis of the dressed DOS by 
examining the appearance of phonon peaks in the sim- 
plest case of vanishing fi and A. As discussed above, 
just as E picks up the oscillations in the bare DOS, 
the renormalized DOS picks up those in E. Figure [2] 
shows the general effects of this renormalization. Here 
we use a fairly weak coupling of A = 50meV, which al- 
lows us to easily compare the renormalized DOS (solid 
black curve in the top frame) to the bare DOS (dashed 
red curve). We also show the self-energy, in the lower 
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FIG. 3. (Color online) The sum of states J Q " N(lo')cIlu' as a 
function of energy for B — 15. 2T (solid black curves) and 
B — (dashed red curves) in both the bare (left) and dressed 
(right) cases. The remaining parameters are A = 250meV, 
uie = 200meV, and F = 5meV. The number of states is writ- 
ten in units of 10 5 meV 2 jV . 



two frames, to make transparent its manifestation in the 
dressed DOS. For small w, -£i(w) = wA cff + 0(w 2 ), 
where A cff = — 7^(0) is the effective mass renormal- 
ization parameter. As seen in the middle frame of the 
figure, this linear behavior remains a good approxima- 
tion in the window (—we, we). And since £2 vanishes in 
that window, the only effect of renormalization therein 
is to draw the Landau levels closer together, due to 
the linear change in Si. For the parameters used here, 
A cff = 0.04 and the n = ±1 LLs, located at solutions to 
M±x = E±x - Ei(P ±1 ) ~ E±x + \ eS E±x, are appropri- 
ately shifted down by a factor of 1 + A eff . This shift from 
M±x to E±x = M±x/(1 + A cff ) can be interpreted as a 
renormalization of the Fermi velocity in M±\. Outside 
the window, above the phonon energy, the levels are still 
drawn closer together but also broadened, as E 2 becomes 
finite with the opening of a new scattering channel. More 
notably, new peaks arise in the DOS, indicated by blue 
arrows and labeled with integers m p . As described in the 
previous section, these peaks lie at M n ± we, the same 
locations as the peaks in E. Specifically, 0+ and l p lie 
at we and M\ +we, respectively, and 0~ and —l p lie at 
—we and M_i — we- 



We can also examine how these features in the 
DOS manifest in the sum of states between and w, 
f£ N(uj')duj' , as shown in Fig. [3] On the left, we see that 
just as in the DOS, the magnetic field causes the bare 
B > sum (solid black curve) to oscillate about the bare 
B = sum (dashed red curve). On the right, we see the 
same general behavior in the dressed sums, except for a 
disruption around the phonon frequency. Hence, in both 
cases, adding a magnetic field does not alter the total 
number of states. 



B. Locations and splitting of dressed Landau levels 

In Fig. g] we show more precisely how and where the 
two sets of peaks in the DOS appear, now choosing a 
stronger, more realistic coupling. We also choose a small 
phonon frequency, such that we can see the first phonon 
peak clearly, with no nearby Landau levels. The upper 
frame shows the DOS and the lower shows a plot of w — 
Ei(w), both in solid black. The energies of the first three 
phonon peaks, labeled 0+ , l p , and 2 p in the upper frame, 
are indicated by vertical dashed red lines. We see in 
the bottom frame that these energies correspond to the 
center point of an oscillation in w — Ei(w). As we can 
see, the first peak is highly asymmetric, since no phonon 
peaks can occur below we- It is the image of half of the 
n = level from the bare DOS; the other half would lie 
at —we- 

However, the plot of w— Ei(w) is primarily useful in de- 
termining the locations of the Landau levels, rather than 
the already-known locations of the phonon peaks. To this 
end, in the lower frame we also plot the bare energies M n , 
shown with horizontal dotted lines; since the renormal- 
ized LLs lie at the solutions to E n — £i(P n ) = M n , their 
renormalized energies E n occur at the intersections of 
these dotted lines with the solid black curve. The ener- 
gies E n are then indicated with dashed blue vertical lines 
passing through the intersection points. Looking at the 
top frame, we see that those intersection points correctly 
determine the locations of peaks in the dressed DOS. We 
also see that two of the LLs, n = 3 and n = 5, are 
split. This occurs due to oscillations in Ej, which allow 
multiple intersections to occur for a given M n . Because 
a phonon peak always occurs at the center of an oscil- 
lation, the LLs will always be split such that one peak 
occurs to the left of a phonon peak and one to the right. 

This splitting occurs only if a Landau level is suffi- 
ciently close to a phonon peak. More precisely, it oc- 
curs only when the bare level at M n lies within the 
amplitude of an oscillation in w — Ex(w) + u (includ- 
ing the possibility of finite \x for generality). Denoting 
the amplitude of the mth oscillation, occurring about 
P m — Ei(P m ) + (J,, by <$Ei(P m ), we can write the con- 
dition for the nth LL to be split about the mth phonon 
peak as \M n - [P m - £ x (P m ) +/*]| < SE x (P m ). We can 



straightforwardly derive a formula for d£j from Eqs. ( 10 ) 
and (11), as follows. Near the energy P m , we can treat 



all but the mth peak as a smooth background and look 



at the function 



AB 



alone. This vanishes at 



ttw (w-p m ) 2 +r 

w = P m and has extrema at w = P m ± T. For suffi- 
ciently small r, we can set the arccos factor that multi- 
plies this equal to tt, so the oscillations have magnitude 
<5Ei = 2WT' independent of m. Noting that the only 
effect of B is to create oscillations about the B = case, 
we have Ei(P m ) = Ef =0 (P m ). After similarly noting 
A* = Mo + Ef =0 (0), we can conclude that in general, the 
nth LL will be split in two about P m when M n satisfies 



\M„ 



[P m -A£i(P m ) + /i<,]| <a£i, 



(12) 
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FIG. 4. (Color online) Locations of peaks in the density 
of states, as determined by to — Ei(w), for A = 250meV, 
lje = 50meV, B = 12. 2T, and T = 4meV. Straight lines indi- 
cate the bare Landau-level energies M n (dotted black lines), 
the energies of renormalized Landau levels E n (dot-dashed 
blue), and the positions P m of phonon peaks (dashed red). 
Blue circles mark the intersection of u) — Ei with M n , which 
determines the energies E n ; since these intersections occur 
twice for M3 and M5, the third and fifth LLs are split. Red 
circles mark the center-points of oscillations in E, which de- 
termine the positions of phonon peaks. Those red circles out- 
lined in blue also satisfy the equations oj — Ei = M n (for n = 3 
and n = 5). 



where ASi(P m ) = Sf=°(P m ) - Sf=°(0) is the change 
in Si away from its value at the Fermi energy, and 
5 Ei = ^ s tne magnitude of the oscillation in Ei 

about its value at P m . With this condition, one can de- 
termine whether a level will be split without any numer- 
ical computation, since all quantities are known analyti- 
cally in terms of elementary functions. [The equation for 
Ef =0 (w) can be found in Carbotte et alW\ In addition, 
the bare chemical potential /io need not be known, since 
its appearance in P m cancels its explicit appearance in 
Eq. pi 



We illustrate the effect of changing T in Fig. [5j which 
shows a small portion of Fig. [4] but now with two values 
of T. The dashed black curves are identical to the curves 
in the previous figure, while the solid green curves are 
for a halved value of T. Halving V doubles <5Ei, which 
results in more Landau levels being split. Here the black 
curve, with T = 4meV, intersects only the n = 3 bare 
LL energy in the lower frame, so only the n = 3 level is 
split in the DOS in the upper frame. The green curve 
intersects the bare energies M 2 through M 5 , leading to 
four split levels in the DOS. Note that the intersection 
with M5 is glancing, such that the n = 5~ level is not 
apparent, instead only influencing the asymmetry of ad- 
jacent peaks. Similarly, the intersections on the side of 



FIG. 5. (Color online) The effect of F on the number of 
peaks in the density of states, showing a small portion of 
Fig. [4] where the dashed black curves here are identical to 
the solid black curves therein, and the solid green curves are 
for a halved value of F. The curve with V — 2meV intersects 
two more M n 's than does the curve with F = 4meV, creating 
two more split Landau levels. 




FIG. 6. (Color online) A schematic depiction of the locations 
of energy levels on the Dirac cone for three increasing values 
of B. Solid colored lines labeled with integers n or rv^ indicate 
Landau levels. Dashed black lines labeled with integers m v 
or indicate energies of phonon peaks. 



the oscillation nearest P m (e.g., the unmarked intersec- 
tions of M2 and M4 with the green curve in the lower 
frame) do not result in a new peak, because those peaks 
disappear under the phonon peak. 

One should note that this analysis does not apply to 
the logarithmic peaks at ±u>e- Referring to Eq. (11), we 
see that except in the case when a peak at P m coincides 



with a logarithmic peak, the logarithmic peak vanishes 
for r — > rather than diverging as the other peaks do. 
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FIG. 7. (Color online) Top: Density of states as a function of 
energy for a sequence of magnetic field strengths, with curves 
offset by 500meV for clarity. Black arrows point to the first 
two sets of phonon peaks, at ±u>e and ±(u>e + Mi). Integers 
label the first five sets of Landau levels. In all cases, uje = 
200meV, A — 250meV, and F = 5meV. Bottom: Locations 
of peaks in the DOS as a function of V~B, with the same 
parameters. The solid curves labeled E n ± show the positions 
of the renormalized Landau levels. The thin dashed lines 
show M n /(l + A ), where Ao = A cff (B = 0). The thick black 
dashed lines show the first two phonon peaks. Inset: the 
fractional change in A cff away from Ao as a function of B, 
for uje = 150meV (dashed red curve), uje = 200meV (solid 
black), and uje = 250meV (double-dot dashed blue). 



same rate as the corresponding LLs, ±m, of which they 
are images. So in particular, the phonon peaks Op remain 
stationary, since they are images of the n = LL, which 
is pinned to the Dirac point. Because the phonon peaks 
arc shifted outward by relative to the bare LLs, the 
LLs pass through them as the magnetic field increases, 
causing the levels to split before reconverging at higher 
B. In the schematic, we see the ±3 levels split about 
the 0^ peaks in the first frame, then reconverge in the 
second; the ±2 levels split about the Op peaks in the sec- 
ond, then reconverge in the third; and the ±5 levels split 
about the ±l p peaks in the third. 

Figure [7] displays these effects more precisely. In the 
upper frame we show a sequence of DOS curves for a 
range of magnetic fields. We see both the Landau lev- 
els and phonon peaks move outward as B increases, with 
the exception of the zeroth level at the Dirac point and 
the Op phonon peaks at plus or minus the phonon fre- 
quency. As the LLs approach the Op peaks, they are 
visibly disrupted, splitting into two peaks of unequal 
width and then merging back into single peaks once they 
are sufficiently far out. In the lower frame we show 
(with solid curves) the energies of the first five (posi- 
tive n) LLs as a function of \J~B. For comparison, we 
also show (with dashed lines colored to correspond to 
the solid ones) the straight lines M n /(1 + A ), where 
Ao = A cff (B = 0). When sufficiently far from any phonon 
peaks, the LLs closely follow the straight lines. Near 
the phonon peaks, they deviate from this behavior and 
split in two, as expected. Note that we could calculate 
curves M„/(l + A cff ), where the renormalization factor 
1/(1 + A cff ) depends on B. This would more closely track 
the true renormalized curves. To determine the magni- 
tude of the correction, in the inset we plot X eS for three 
values of cue- As we would expect, bringing the phonon 
closer to the Fermi energy increases its effect on the slope 
of Sr there, increasing the change in A with B. How- 
ever, for the value uje — 200meV used in the remainder 
of this paper, all changes are less than 5% for a very large 
range of field strengths. So, except near a phonon peak, 
M n /(l + \a) is an accurate approximation for the energy 
of the nth renormalized LL. This result, in combination 
with the energies of the phonon peaks given in Table [I] 
and the splitting condition in Eq. ( fl2"| ), allows us to cal- 
culate the approximate positions of all peaks in the DOS 
in terms of elementary functions. 



The foregoing results have all been for a fixed mag- 
netic field. We now consider the effect of varying that 
field. The general trend is shown schematically in Fig. |6j 
which displays the position of the peaks on the Dirac 
cone for three increasing values of B. As the magnetic 
field is increased, the dressed Landau levels (shown as 
colored lines labeled with integers n or n ± ) are pushed 
further away from the Dirac point. The phonon peaks 
±TOp (shown as dashed black lines) move outward at the 



C. Broadened phonon spectrum 

Thus far, our results have been restricted to the sim- 
ple Einstein model. In a real physical system, the phonon 
distribution will be broadened, and this broadening can 
significantly affect the features of the DOS around the 
phonon energy. As an example, we conside r the trun- 
cated Lorentzian distribution of Dogan et al.j 44 * 45 * given 
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FIG. 8. (Color online) Density of states and self-energy in 
the Einstein model (dashed red curves) and in the truncated 
Lorentzian model (solid black curves) with a magnetic field 
B — 15. 2T. The Landau levels that are split about oje in the 
Einstein model are flattened by the broadening of the phonon 
in the Lorentzian model. In the Einstein model we use A = 
250meV and loe = 200meV; in the Lorentzian, luq — 200meV, 
A' = 555meV, 5 = 15meV, and S c = 30meV. The value of A' 
is chosen to ensure that A eff = 0.19 in both models. In both 
cases, F = 5meV. 



by 

a " F ^ = - \l Wa^ - actI W c -\v-w \). 

7T \[y — woj + o£ + o z \ 

(13) 

This is a Lorentzian of width S centered at loq and trun- 
cated at uo ± S c . In Fig. [8] we present results for this 
model together with the Einstein model. To compare the 
results, we set ujq = <^e and choose A' such that A off is 
identical for both models. We see that in the Lorentzian 
model, the magnitude of peaks in the self-energy are re- 
duced from that of the Einstein model. More signifi- 
cantly, the peak at lje is smoothed, losing the distinctive 
one-sided character it has in the Einstein model. This 
smoothing alters the peak structure around loe in the 
DOS, merging the split Landau level with the phonon 
peak to create a flattened structure. As discussed in our 
companion paperj^ this flattened structure can be seen 
in experimental data. 



D. Multiple iterations 

In addition to restricting ourselves to an Einstein 
model, we have thus far restricted ourselves to one it- 
eration of our self-consistent system of equations. In 
Fig. [9] we present the DOS and self-energy after iterat- 



FIG. 9. (Color online) Density of states and self-energy after 
one iteration (dashed red curves) and after iterating to con- 
vergence (solid black curves), with parameters B — 15. 2T, 
A = 250meV, lo e = 200meV, and F = 5meV. 



ing to convergence. The slightly more complicated set 
of peaks that appears in the DOS can still be divided 
into Landau levels, which solve E n — Hi(E n ) + /i = M n , 
and phonon peaks, at which £2 has its maxima. We see 
that the converged function £2(w) contains a perfect im- 
age of the dressed, rather than the bare, DOS, such that 
phonon peaks now lie at E n + uje rather than M n + uje- 
This result follows obviously from Eq. [6] Perhaps less 
obviously, phonon peaks representing multiphonon pro- 
cesses also arise, at energies E n + jcos for integers j > 2. 
These peaks occur because the factors N(u> ± u>e) in 
the self-energy pick up not only the LLs, but also the 
phonon peaks in the preceding iteration. So, for exam- 
ple, X4 2) (w) °c N^{cj - lu e ) will have a phonon peak 
at 2loe, corresponding to the phonon peak in iV*- 1 ' at 
loe- And each further iteration produces a new set of 
phonon peaks shifted further away from the Fermi en- 
ergy by wg. These new phonon peaks can cause further 
splitting of LLs. Moreover, the phonon peaks that were 
already present in the first iteration (shown in red) are 
shifted down in energy, causing split LLs at lower en- 
ergies. Together, these effects of iteration cause the fea- 
tures in the DOS to become progressively more indistinct 
at energies > E\ +uje- 

However, at energies below this, which are more rele- 
vant to current experiments, the only effect of iteration is 
a slight increase in A; the qualitative features are wholly 
unaffected. And for smaller values of A, the effects are 
further reduced. For example, for the value A — lOOmeV 
used in our fit to experimental data^, iterating to con- 
vergence alters the parameters obtained from the fit by 
less than 1%. Given this fact, in the following sections 
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we will continue to use the Einstein model and only one 
iteration, since that affords the cleanest description of 
features in the DOS. But one should keep in mind two es- 
sential modifications that would occur in a more complete 
treatment: (i) near phonon peaks, a broadened phonon 
spectrum will smear features in the DOS; and (ii) in a 
converged solution, phonon peaks are displaced by 
relative to the renormalized Landau levels rather than 
the bare ones. 



IV. DENSITY OF STATES FOR FINITE 
CHEMICAL POTENTIAL 

A. Jumping phonon peaks 

When the chemical potential is finite, many new fea- 
tures appear in the DOS. First among these is a jump 
in the location of the phonon peaks as B is varied. Re- 
ferring back to Table |TJ we see that whenever a phonon 
peak at P m = M m — /iQ ± uje would fall in the window 
(— uj e , we), its position is altered to P m = M m — hoT^e, 
causing it to jump from one side of the window to the 
other. This discontinuity in the positions of peaks is seen 



in the bottom frames of Fig. 10 which show plots of P m 
as a function of \ B for the two cases fio > and « < 0. 
Here we see that as B is varied, the phonon peaks jump 
from one side of ±we (dotted black horizontal lines) to 
the other side of T^e. More precisely, in the left plot, 
when B > //q/2 none of the curves have reached ±uj e , 
and they are given by — « (black), M±\ — [i ±uje 

(red), M±2 — Ho±^ E (green), and M ±3 ~ fi Q ±uj E (blue), 
where upper signs correspond to upper curves and lower 

2 

to lower. When B reaches the nth energy above uj e , 
given by M n — /iq+uje, reaches uje and then jumps discon- 
tinuously to below — u>e, continuing on as M n — Ho — u>e 
at lower B. The energies in the right plot differ from 
those in the left only by a change of sign. 

This phenomenon allows us to single ou t three cases, 
shown in the top three frames of Fig 
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-V2B 



« > V2B, 
In the first 
and at least one 
> occur to the left of the window 



-V2-B < « < v2-B, and fio < 
case, all phonon peaks with m p < 
phonon peak with m p 
(~loe,uje), while the remaining phonon peaks lie to the 
right; for example, in the top frame of the figure, l p lies 
to the left. In the second case, all phonon peaks with 
trip < lie to the left and all with m p > lie to the 
right. In the last case, at least one phonon peak with 
Tftp < lies to the right; for example, — l p in the third 
frame. The position of the zeroth level relative to the 
window depends only on the sign of Ho'- it lies to the 
left for hq > and to the right for Ho < 0. As seen in 
the preceding section, for fio = 0, the peak is split into 
two, with one at — uje an d one at w^. In all cases, the 
logarithmic peaks lie at ±u>e- 

Figure [TT] shows this jumping behavior schematically 
on the Dirac cone. Just as in the corresponding schematic 
for [i — 0, shown in Fig. [6j solid colored lines labeled 
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FIG. 10. (Color online) Top: Density of states as a func- 
tion of energy for three finite values of /i: 150meV (top), 
50meV (middle), and -150meV (bottom), chosen to lie in the 
intervals fio > v2B, —y2B < fj,o < v2B, and fio < —V2B, 
respectively. Arrows point to the locations of the phonon 
peaks. In all cases, B = 15. 2T, uj e = 200meV, A = 250meV, 
and r = 5meV. Bottom: the locations of the phonon peaks 
for no — 150meV (left) and fio = — 150meV (right). The col- 
ors of the curves match those of the arrows in the top set of 
plots. Dotted lines indicate ±ub. 



with integers n are Landau levels and dashed black lines 
labeled with integers m p are phonon peaks. The cones 
are now shifted by —fj, relative to the Fermi energy (the 
filled region below the Fermi energy is shaded), so in 
the upper set of cones, for /x > 0, the cones are shifted 
downward, and in the lower set, for u < 0, the cones 
are shifted upward. As in the earlier schematic, both the 
LLs and the phonon peaks move outward from the Dirac 
point as B increases, and the LLs split when passing 
through a phonon peak. However, now as we increase B, 
the locations of the phonon peaks vary discontinuously. 
On the upper set of cones, we see first 2 p and then l p 




Bi B 2 > Bi B 3 > B 2 



FIG. 11. (Color online) The arrangement of energy levels on 
the Dirac cone in the cases u, > (top) and /i < (bottom) for 
three increasing values of magnetic field. The shaded region 
indicates filling up to the Fermi energy. Landau levels are 
indicated by solid lines labeled with integers n or n , and 
phonon peaks by dashed lines labeled with integers m v or 
nip. 

jump from below — uje to above uie- The lower set is the 
mirror of the upper, so as we increase B, we see first — 2 p 
and then — l p jump from above uie to below —uie- 



B. Large chemical potential and split zeroth level 

When > lue, the n = Landau level lies outside 
the window (—lue, ^e), introducing several new features. 
This scenario is shown in Fig.[l2j where we show both the 
bare and dressed DOS. We also show the B = curves 
for comparison. For B = 0, the DOS at the Dirac point 
(at — ^io in the bare case and slightly to the left of — // in 
the dressed case) is pronouncedly lifted by renormaliza- 
tion. In the finite- B case, this occurs as well, and we see 
here that the valleys around —fi are no longer global min- 
ima, as they were in the bare case. But additional, even 
more marked changes occur in the finite- B case: there is 
no single peak near — \x that indicates the Dirac point. 
Recall that since Si is effectively unchanged by the pres- 
ence of a magnetic field in the window (— uie, ^e), the 
chemical potential [i = + Si(0) is also unaffected by 
the field, so — fio and — // refer to both the B = and 
B > cases. 
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FIG. 12. (Color online) Bare (thin dashed blue) and renor- 
malized (solid black) density of states and the real part of 
the self-energy in the case of [i > u>e- Arrows point to the 
(negative of the) bare and renormalized chemical potentials 
fi — 340meV and fi = 291meV, respectively. The remaining 
parameters are uie = 200meV, A — 250meV, B — 15. 2T, and 
r = 5meV. We also display the bare (dot-dashed magenta) 
and renormalized (thick dashed red) curves in the B = case. 



We examine this situation more closely in Fig. [T3j 
which shows the region around —fi in more detail, label- 
ing all the peaks in the region. We see that the three-peak 
structure near —\i consists of a phonon peak, 3 P , and a 
split zeroth LL. Hence, when fi lies outside the window 
(— uje,uje), the zeroth LL loses much of its significance: 
it may be split just as any other LL is when sufficiently 
close to a phonon peak, and it is not pinned to the Dirac 
point. (Also note the distinct character of the peak at 
— uje, arising from the logarithmic terms in Si, which 
for these particular parameters is clearly visible because 
none of the phonon peaks' energies P m lie near —uie-) 

Since there is no single zeroth LL, we have no clear way 
of identifying a Dirac point here. In the B = case, we 
can define the Dirac point even when we include interac- 
tions. The bare band energies — ±VFk— no become the 
dressed energies Ek, given by Ek — ^i(Ek) — ±vpk — /i. 
We can then define the Dirac point as the point of 
k = on these curves, such that its energy ui^ satis- 
fies LUd — Si(wd) = — fJ" Although the dispersion relation 
is only approximately linear in this case, it is sufficiently 
close to linearity for the Dirac point to remain meaning- 
ful. In the finite- B case, however, the linear bands vanish, 
replaced by discrete LLs. Even when neglecting interac- 
tions, this seems to make the definition of the Dirac point 
ambiguous. We illustrate this situation in Fig. [TIJ which 
shows color maps of the spectral density as a function of 
k and u>; the black lines are intended only to guide the 
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FIG. 13. (Color online) A small portion of Fig. 12 focusing 
on the region around — fi. Here we see that the three peaks 
near — fi consist of a phonon peak and the split n = Landau 
level. None of the three lies at — /x. 
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eye. The spectral density here is given by^l 

oo 

tr[ 7 °A(fc,u;)] = 2e-^ fc2 / B £(-l) n 
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X — 
7T 



L n (2v 2 F k 2 /B) - L n ^(2v 2 F k 2 /B)] 

r-s 2 



r-s 2 



(w - Si + n + M n y + (r - e 2 ) 2 



(14) 



where the spectral function A(k,cd) is a 4-by-4 block- 
diagonal matrix with one block for the K point and one 
for the K' point, 7 = (73(g) 0-3 is the Dirac matrix (where 
(T3 is the Pauli matrix), k = \k\, and h = 1. As we see in 
the upper frame of the figure, in the bare case the LLs 
are arranged on two fairly definite cones that intersect at 
the n = level. So although there are no longer linear 
bands, we can still define a Dirac point based on this 
conic arrangement of levels, and we can identify the n = 
level with the Dirac point. But in the dressed case, shown 
in the lower frame, the conic arrangement is radically 
distorted, particularly by the splitting of the n = level. 
This leaves considerable ambiguity in defining the cones 
or a Dirac point. However, we can see from the figure 
that two intersecting lines alone cannot be made to pass 
through the peaks in the spectral function. In that sense, 
we can say that the splitting of the zeroth LL corresponds 
to a splitting of the Dirac point into two. 



FIG. 14. (Color online) Color maps of the bare (top) and 
dressed (bottom) spectral density tr[7°A(fc, co)] in the (k,Lo) 
plane for no = -340meV, uj e = 200meV, A = 250meV, B = 
15. 2T, and F = 3meV. The plane is extended to negative 
values of k in order to show the cone-like structure. 



Although this picture of the spectral density would be 
altered if electron-electron interactions were taken into 
account as well, we expect it to remain qualitatively valid. 
In the figure we see each energy level broken into a series 
of islands of strong spectral weight, due to the Laguerre 
polynomials in the spectral function. (These islands arise 
from the multiplication of the two sublattice compo- 
nents in the wavefunction in a magnetic field, the K- 
point piece of which reads^-3 ^ X^n(V|n|' ^M+i) m the 
(A, B) sublattice basis, where ip n is the harmonic oscil- 
lator eigenfunction.) Below —uje, we see that the islands 
are smeared at the energies of phonon peaks; above uje, 
the islands become indistinguishable because their spac- 
ing is on the same scale as the phonon peaks. If electron- 
electron interactions were accounted for, the particle-hole 
excitations would also be localized into islands,^ with a 
spacing in w equal to the difference between two LLs and 
a spacing in k determined by the wavefunction of the 
lower level, which would cause smearing at those loca- 
tions. But unlike in the zero-field case, the excitations 
are gapped out at the Dirac point, so the splitting of the 
n = level due to electron-phonon coupling could be 
distinguished. 
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FIG. 15. (Color online) Density of states with a finite en- 
ergy gap A = lOmeV for three values of fio'- 150meV (top), 
OmeV (middle), and -150meV (bottom), chosen to lie in the 
intervals /Uo > A, —A < fio < A, and fio < — A, respectively. 
Red arrows point to the locations of the phonon peaks around 
±u>e, shown in the insets. For —A < fio < A, these peaks lie 
at ±(lje + A) — fio; for fio > A, at — uje±A — fio; for fio < —A, 
at uje ± A — fio. In all cases, uje = 200meV, A = 250meV, 
B = 15. 2T, and T = 5meV. 



V. OTHER EFFECTS CREATING NEW PEAKS 

In this section, we analyse two more ways in which 
phonons can induce new peaks in the DOS. First, we 
discuss how coupling to a phonon causes new structures 
to arise from the splitting of the zeroth LL by an energy 
gap. Next, we discuss the effects of an acoustic phonon 
modeled by a Debye distribution. 



A. Finite energy gap 

Recent experimental data has suggested the existence 
of an energy gap A at the Dirac point in grapheneP The 
cause of this gap is unknown, but may be related to sub- 
lattice asymmetry! 32 ! 33 ! As we show in this section, the 
signature of a gap appears not only in splitting the n = 
Landau level, but in the images Op of that split level in 
the phonon peaks at ±cj£. From Table [IJ we see that the 
locations of the Op levels depend on the relative size of 
A and fiQ. The three possible cases are shown in Fig. 15 
with the locations of the peaks indicated by red arrows. 
If /i lies between —A and +A, as shown in the mid- 
dle frame of the figure, then the image of the bare level 
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FIG. 16. (Color online) A density of states curve displaying 
a Landau level (n = —3) split in three by a phonon peak that 
is itself split in two by an energy gap. Here fio = 15meV, 
A = lOmeV, A = 250meV, u) E = 200meV, V = 3meV, and 
B = 15. 2T. 



n = + at oj = A — fi appears to the right of uje, as the 
0+ level ai uj = A — fio+UE, while the image of the bare 
level n = 0~ at u) = —A — fiQ appears to the left of —uje, 
as the 0~ level at oj = — A — fiQ — cue- If Mo lies outside 
this range, then both the 0+ and 0~ levels appear on the 
same side of the window to the right of uje, 

at ±A — /io + t-^E, if Mo > A; or to the left of —loe, at 
±A — )1q — uje, if Mo < —A. These two cases are shown 



in the top and bottom frames of Fig. 15 respectively. 



If the Op peaks both lie on the same side of the window 
(— loe,uje), then any sufficiently close Landau level will 
be split into three peaks rather than two. We illustrate 
this in Fig. [T6j where the LL n = —3 is split in three. 
Because the bare n = ± levels have equal magnitude, 
their images Op do as well, creating the two neighbouring 
oscillations of equal magnitude seen in ui — Si + /i. For 
sufficiently small A, this means that the splitting of a LL 
into three is generic, since a line of constant M n will have 
to intersect both oscillations if it intersects cither. 



B. Acoustic phonon 

Recent STS experiments have also observed shoulders 
and distinct secondary peaks on the Landau levelsP^l 
While some or all of these peaks may be caused by the 
STM tip deforming the graphene sheet they could po- 
tentially be caused by coupling to an acoustic phonon 
or some other low-energy boson. To consider this case, 
rather than the Einstein model, we adopt a Debye distri- 
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FIG. 17. (Color online) The bare density of states (dashed 
red curves) and the density of states renormalized by an 
acoustic phonon (solid black curves), for three values of fio' 
15meV (top), OmeV (middle), and -15meV(bottom). In all 
cases, B = 15. 2T, F = 5meV, and the parameters of the De- 
bye distribution are ojd = 20meV and A' — 0.03125meV~ 2 . 



bution for the electron-phonon spectral density, 

a 2 F{v) = A'is 2 6(uj d - v). (15) 



Figure [17| shows the DOS for three cases: Ho > 0, = 0, 
and no < 0. When the chemical potential vanishes, the 
acoustic phonon manifests as a hump on either side of 
the zeroth LL, along with humps at M n ±u>d- When the 
chemical potential is finite, only one hump appears near 
the zeroth LL, to the left for /x > and to the right 
for Ho < 0, in the same manner that Pq appeared to the 
right or left in the case of the Einstein model at finite //. 
The other humps are unaffected, only shifted by — ho to 

From this, we see that coupling to a Debye distribution 
has qualitatively the same effect as coupling to a single 
Einstein mode, introducing a new set of peaks that is 
shifted relative to the Landau levels by the characteristic 
phonon frequency (in this case, the Debye frequency). 
However, since the Debye distribution is finite at all 
nonzero frequencies below wr>, in the Debye model the 
interval (-Wd,Wd) is 'softer' than the analogous interval 
(— loe, oje) in the Einstein model. By that we mean that 
phonon peaks can appear within the window, but with 
greatly reduced magnitudes, and LLs within the window 
will be more affected by the phonon. But the Debye fre- 
quency is likely to be small compared to the scale of the 



LLs in experiment, so the window will contain few levels 
in any case. 



VI. CONCLUSION 

In this paper, we have shown how in the presence of 
a magnetic field B, electron-phonon coupling induces a 
new set of peaks in the electronic density of states in 
graphene, over and above the Landau levels in the bare 
DOS; furthermore, we have shown how these new peaks 
affect the locations of the renormalized LLs themselves. 
We can understand the location of the new peaks, which 
we have called phonon peaks, by noting that the electron- 
phonon self-energy £ is associated with an electron (or 
hole) scattering between states of energy u> and u> — uie 
(or uj + oje), where u>e is the phonon energy. This means 
that the self-energy at oj contains a direct imprint of the 
DOS at cj ± uje, leading to peaks in £ at energies dis- 
placed by iuiE from the LLs. Renormalization then in- 
troduces these phonon peaks into the dressed DOS. Since 
this scattering process is possible only for \oj\ > uje, the 
phonon peaks always lie outside the window (—uje,<-^e)- 
The presence of these new peaks alters the behavior of 
the LLs as well. When a LL lies far from a phonon peak, 
renormalization simply scales its energy by a factor of 
1/(1 + A ) (equivalent to a renormalization of the Fermi 
velocity) , where A off is the effective mass renormalization 
parameter. However, the nearer a LL lies to a phonon 
peak, the more it deviates from this behavior; and if it 
is sufficiently near, it splits in two about the phonon 
peak. Because the phonon peaks lie outside the win- 
dow (— uie, we), LLs within that window are left almost 
entirely intact, simply shifted by the renormalization of 
the Fermi velocity. 

This picture is made more complex in the case of a 
finite chemical potential, when the Landau levels and 
phonon peaks are no longer centered on the Fermi en- 
ergy. Because no phonon peaks can occur in the win- 
dow (—uje,i^e), varying B or fi continuously causes the 
phonon peaks to jump discontinously from one side of 
the window to the other. And if the chemical potential 
is greater than the phonon frequency, then the zeroth 
LL lies outside the window, falling amongst the phonon 
peaks. This means that it can be split in two, losing its 
signifcance as an identifier of the Dirac point. 

We have also discussed how new peaks arise in the 
presence of an energy gap at the Dirac point or a cou- 
pling to an acoustic phonon. In the first case, the energy 
gap splits the bare zeroth LL, leading to a split zeroth 
phonon peak P ± . Depending on the value of the chem- 
ical potential, the peaks Pq± can both lie to one side of 
the window (— lue^e), or one can lie to either side of 
it. If both lie on the same side, they can split a nearby 
LL in three. In the second case, coupling to an acoustic 
phonon modeled by the Debye distribution has the same 
qualitative effect as coupling to an Einstein mode, intro- 
ducing phonon peaks at M n ± uje> . In the case of charge 
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neutrality, this leads to two phonon peaks centered about 
the zeroth LL; for finite chemical potential, there is only 
one zeroth phonon peak, lying to the right or left of the 
zeroth LL. 

Together, these analyses reveal a rich electronic struc- 
ture that emerges from an interplay between Landau level 
quantization and electron-phonon coupling in graphene. 
Because the features in this structure are tied to a fixed 
frequency that doesn't depend on magnetic field, they 



should be observable in experiment even in the presence 
of complicating factors such as electron-electron interac- 
tions, finite size effects, and substrate effects. Indeed, 
some of the features have already been observed in STS 
dataj^ others, such as the hopping of phonon peaks at 
finite /i, should be readily observable, while precise ob- 
servations of some, such as the behavior of split Landau 
levels, may require higher quality samples. And future 
studies of different quantities, such as optical conductiv- 
ity, should yield similarly rich results. 
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